Renormalization-group study of gate charge effects in Josephson-junction chains 
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We study the quantum phase transition in a chain of superconducting grains, coupled by Joseph- 
son junctions, with emphasis on the effects of gate charges induced on the grains. At zero temper- 
ature the system is mapped onto a two-dimensional classical Coulomb gas, where the gate charge 
plays the role of an imaginary electric field. Such a field is found relevant in the renormalization- 
group transformation and to change the nature of the superconductor-insulator transition present 
in the system, tending to suppress quantum fluctuations and helping establish superconductivity. 
On the basis of this observation, we propose the zero-temperature phase diagram on the plane of 
the gate charge and the energy ratio. 
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In arrays of ultra-small superconducting grains, the 
charging energy can be dominant over the Josephson cou- 
pling energy^ giving rise to crucial effects of quantum 
fluctuations.!!] In particular recent advances in fabrica- 
tion techniques make it possible to control the physical 
parameters of the arrays, providing a convenient model 
system for the study of quantum phase transitions 2 ] be- 
tween superconducting and insulating phases. In such an 
array of ultra-small junctions, frustration can be intro- 
duced not only by applying magnetic fields but also by 
inducing gate charges; these control the number of vor- 
tices and charges (Cooper pairs), leading to interesting 
phasa-,tpansitions and dynamic responses in two dimen- 
sions.Hti In one dimension, where quantum fluctuations 
should be more important, the quantum phase transition 
has been studied in the absence of gate charges.B In such 
a one-dimensional (ID) array, i.e., a chain of Josephson 
junctions, gate charges on grains, breaking the particle- 
hole symmetry, cart-generate the persistent voltage in the 
appropriate regimccl and affect the quantum phase tran- 
sition substantially, which has beea studied analytically 
via the mean-field, approximation^] and via the pertur- 
bation expansion!^ as well as numerically via quantum 
Monte Carlo simulations.!!!! In the analytical approaches, 
however, the second gives accurate results only in the 
weak-coupling limit, i.e., for sufficiently small values of 
the Josephson coupling energy, whereas the first is not ex- 
pected to be reliable in one dimension, where fluctuations 
arc too strong to be neglected. Both approaches fail to 
capture the essential physics rfitbe Kosterlitz-Thouless- 
Berezinskii (KTB) transition)! 2 ]!! 3 ] present in the system 
without gate charges. They not only disallow one to ex- 
plore the detailed nature of the transition but also yield 
phase boundaries near the particle-hole symmetry line in. 
rather large discrepancy with those from simulations .Ell 
It is thus desirable to study the gate charge effects in 
an analytical way beyond the above approaches, which 
discloses the nature of the transition. 

This paper presents the renormalization-group (RG) 



study of the quantum phase transition in a chain of 
Josephson junctions with gate charges. At zero tempera- 
ture, the system is mapped onto a two-dimensional (2D) 
classical Coulomb gas, where the gate charge plays the 
role of an imaginary electric field. This in general leads 
to the superconductor-insulator transition, driven by vor- 
tices, as the ratio of the Josephson coupling energy to the 
charging energy is varied. Here the field associated with 
the gate charge is found to grow under the RG trans- 
formation and tend to suppress quantum fluctuations, 
helping to establish superconductivity in the system. On 
the basis of this observation, the zero-temperature phase 
diagram is proposed on the plane of the gate charge and 
the energy ratio. 

We consider a ID array of N superconducting grains, 
each coupled to its two nearest-neighboring grains via 
Josephson junctions of strength Ej. Each grain is char- 
acterized by the self-capacitance C, leading to the charg- 
ing energy scale Ec = 4e 2 /C, and a gate charge Q can 
be induced on it externally, e.g., by applying gate volt- 
age V with respect to the ground, giving Q = CV . Such 
a Josephson-junction chain is described by the Hamilto- 
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where the lattice constant has been set equal to unity, 
the energy has been rescaled in units of the Josephson 
plasma frequency huj p = t/EcEj, the square root of the 
energy ratio K = J Ej /Ec corresponds to the effective 
coupling, and q = Q/2e is the uniform gate charge in 
units of the Cooper pair charge 2e. The number n x of 
Cooper pairs and the phase <p x of the superconducting 
order parameter at site x are quantum-mechanically con- 
jugate variables: [n x , 4> x '\ = i5 xx > . Due to the periodicity 
and symmetry in q, we need to consider only the range 
0<q<l/2. 

Following the standard procedure,0H we introduce the 
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imaginary time r running on the interval [0,/3], and 
write the partition function in the imaginary-time path- 
integral representation: 
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with the action 
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where the inverse temperature has been rescaled accord- 
ing to h(VpP — > f3, and the (imaginary) time slice 5t has 
been chosen to be unity (in units of hu p ). Here d x and 
d T denote the difference operator with respect to the 
position x and to the imaginary time r, respectively: 
d x 4>x,r = 4>x+i,t ~ 4>x,t and d T (j) XiT = (j) x ,r+l - 4> x ,t, 
and the zero-temperature limit (3 — > oo as well as the 
thermodynamic limit N — > oo is to be taken. We then 
apply the Villain approximation to integrate out the 
phase variables {4> x ,t} and perform the duality trans- 
formationliil to write the partition function in the form 
Z = X){s R } ex P{ — ^[ s ]} with the action 



A[s 



1 

2K 



(4) 



R 



where the integer field sr is defined on the dual lattice 
site R = (X,Y), dual to the original 2D space-time lat- 
tice r = (x,t). The Poisson summation formula allows 
us to decompose sr into a real-valued field Or and a 
new integer field tor describing the vorticity. Shifting 
the field variable in such a way that Or — qX = Or, we 
finally obtain 
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where the tilde sign has been dropped for simplicity. This 
leads, apart from the spin-wave part, to the 2D system 
of classical vortices described by the action 



with the vorticity neutrality condition X^rTOr = 0> 
where G'(R— R') = 2tt[G(0) - G(R-R')] describes the 
lattice Coulomb Green's function with the singular diago- 
nal piece G(0) subtracted and y corresponds to the vortex 
fugacity. Note that the gate charge q plays the role of an 
imaginary electric field acting on vortices. Accordingly, 
the effects of the gate charge on the quantum phase tran- 
sition in the Josephson-junction chain reduces to those of 
the (imaginary) electric field in the 2D Coulomb gas. 

To investigate the quantum phase transition displayed 
by the action in Eq. (^), we now consider the correla- 
tion function, with attention to the effects of the gate 
charge. To the lowest order in the fugacity y, the vortex 
correlation function reads 
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which, like the Hamiltonian in Eq. (|lj), depends periodi- 
cally on q. We also have the average vorticity 

(m ) = -2iy 2 ^2 R ~ 27rK sin 2 ^ x = °> ( 8 ) 

R 

upon summing over the orientation; the gate charge does 
not alter the average vorticity. Renormalization of the 
phase correlation function by such vortices then leads 
to the scaling equations describing the criticaL-behavior. 
The phase correlation function takes the formlla 

S (r-r') ee ( e ^--^)) = g sw ( r -r')g v (r-r') (9) 

with the spin wave part 
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and the vortex contribution given byEf] 

bxg v {r-r') = -jy 2 2 * K cos27rqx [(r • V R > Ro 
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Here Eqs. (Q) and (||) have been used,_and the function 
Mr defined on the dual lattice satisfiestj 



dxu R = -dy[G'(r-R) - G'(r'-R)] 
Oyur = dx[G'(r-R) - G'(r'-R)]. 



(12) 



Replacing the lattice summation over ro by the integra- 
tion in the polar coordinate (ro,#o) and noting [(ro • 
Vr )ur ] 2 = r 2 [cos 2 9 (dx o u Ro ) 2 + sin 2 6> (<9y o w Ro ) 2 + 
2cos9 sm8o(dx UR )(dY UR )} 1 we obtain 
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where Jq and Ji are Bessel functions and we have usedtS This leads to 
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Equation (Q) thus reduces to 
lng(r) = mg> st( ,(r) + ln^(r) 

1 ,G'(r)-2^ 2 y 2 G'(r) / dr r^ 27rif J (2nqr ) 
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which defines the effective coupling K e ff. It is then 
straightforward to derive the scaling equations 
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the third of which demonstrates that the gate charge is 
relevant, growing under the RG transformation: q(£) = 
qe l with q = q(£=0) being the initial (physical) value 
of the gate charge. This allows us to write the scaling 
equations in the form 
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which reduce, for q j=.£Lto the standard scaling equations 
for the XY model.lialii. Accordingly, the system under- 
goes a KTB transitionEj"EJ from the insulating phase to 
the superconducting one as K is increased.!] 

When the gate charge is present (g ^ 0), on the other 
hand, the Bessel function in Eq. (|l6|) first decreases with 
t and can even become negative as the renormalization 
proceeds. The effective temperature K~ l then increases 
less rapidly and may even decrease under RG transfor- 
mation; in this way the gate charge tends to suppress 
quantum fluctuations, helping to establish superconduc- 
tivity in the system. Note also that the periodicity in q 
has been lost in Eq. (|l5|), which results from the contin- 
uum approximation, i.e., the replacement of the lattice 
summation over ro by the integration over vq and 6q. To 
be more accurate, we may avoid the continuum approxi- 
mation and employ the Poisson summation formula: 

^/(r )= J2 J dr ro J de f(r ,e Q ) 

ro — oo 

x exp[27rz7"o(j cos^o + fcsin0o)]< (17) 



which indeed manifests the periodicity, in place of the 
first one in Eq. (|l6|); the second one remains unchanged. 
For q = 0, the terms of j ^ as well as those of k ^ 
give only small quantitative corrections, and we may use 
the anpxnximate scaling equations given by Eq. (|l^) or 
(|16|).EjO On the other hand, for q ^ 0, the contributions 
of nonzero j terms can have crucial effects, and should 
be taken into account. 

We have thus computed numerically the RG flow from 
the above scaling equations for various (initial) values of q 
as well as K, and show the typical RG flow in Fig. |l|. The 
dotted line represents the initial configuration, i.e., the 
physical value of the fugacity, chosen to be y = e - ^ / 2 ) A . 
It is obvious that the Gaussian fixed line (y — 0) is sta- 
ble for K > 2 /it and unstable for K < 2/tt. When the 
initial coupling is strong [K ^> 1), the system thus ap- 
proaches the Gaussian fixed line with finite (renormal- 
ized) values of K, yielding superconductivity; at weak 
(initial) couplings (K < 1), the system flows to the high 
(effective) temperature limit (K^ 1 — * oo) and becomes 
insulating. For intermediate couplings, the non-zero gate 
charge, growing under the RG transformation, can make 
K increase and alter the flow qualitatively: The RG flow 
starting near K ss 2/-7T can eventually turn around, ap- 
proaching the fixed line, as illustrated in Fig. [j]. Note that 
the projections of the flow trajectories onto the K~ 1 —y 
plane are depicted here. In fact the trajectories exist 
in the three-dimensional {K~ l ,y,q) space, and do not 
cross with each other. The minimum initial value of the 
coupling, which yields such turning flow, then gives the 
critical coupling strength K Cl below which the system is 
driven by quantum fluctuations to the insulating state. 
Via extensive numerical computation, we have obtained 
K c for various values of the gate charge q and found that 
K c , starting from 0.748 at q = 0, decreases monotoni- 
cally with q; in this way the gate charge tends to suppress 
quantum fluctuations, helping to establish superconduc- 
tivity in the system. Of particular interest is the case of 
the maximal gate charge (q = 1/2), where the presence 
of degeneracy is believed to-.iestore superconductivity at 
arbitrarily weak couplings Jj'E£l The RG approach, with 
the contributions of several nonzero j terms in Eq. Jl8| ) 
taken into account, indeed indicates that K c reaches its 
minimum at q — 1/2. Unfortunately, however, it is not 
allowed to confirm that K c actually vanishes in this case: 
For large gate charges (q ~ 1/2), the scaling equations 
yield the RG flow going through the region of large values 
of the fugacity (y > 1), for which the scaling equations 
themselves are not reliable. 

The critical coupling strength K c , obtained numeri- 
cally from the scaling equations as above, leads to the 
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phase diagram in Fig. |2j. Note that due to the employed 
Villain approximation, the coupling K in the scaling 
equations cannot be identified directly with y/Ej/Ec 
in the original Josephson-j unction chain system. Com- 
parison of the interaction form in the Villain action and 
the original one shows that the critical coupling strength 
K c = 0.748 for q = corresponds to the critical value 
(Ej/Ec) c ~ 1.035, which sets the overall scale. Accord- 
ingly, for each value of q, the obtained value of K c is 
transformed into the corresponding value of (Ej / Ec) c , 
which gives the boundary between the superconducting 
phase and the insulating one in Fig. |^; manifested is the 
suppression of quantum fluctuations by the gate charge. 
At q = the transition between the two phases is of the 
standard KTB type, driven by vortices in the 2D space- 
time. On the other hand, the nonzero gate charge, which 
breaks the particle-hole symmetry, makes a relevant per- 
turbation and changes the transition in a significant way. 
The critical RG flow corresponding to K c approaches a 
fixed point on the Gaussian fixed line with an arbitrar- 
ily large renormalized value instead of the standard KTB 
value 2/-7T, thus suggesting a transition of the Gaussian 
nature. For large q, the RG flow goes through the re- 
gion of large values of the fugacity, for which the scaling 
equations are not reliable. The corresponding boundary, 
which is thus somewhat speculative, is drawaas a dashed 
line. Note that the perturbation expansion)!!] valid in the 
limit Ej/Ec — > 0, gives accurate results around q « 1/2. 
In this sense the RG approach here complements the per- 
turbation expansion. Here it is pleasing that Fig. |^ agrees 
well with the results of quantum Monte Carlo simula- 
tionsjlil particularly near the symmetry line, where the 
perturbation expansion as well as the mean-field approx- 
imation gives large discrepancy. 

In summary, we have employed the renormalization- 
group transformation method to investigate the quan- 
tum phase transition in a chain of Josephson junctions, 
with attention to the effects of gate charges. At zero 
temperature the system has been mapped onto a two- 
dimensional classical Coulomb gas, which undergoes a 
superconductor-insulator transition as the energy ratio is 
varied. In the absence of the gate charge the transition, 
driven by vortices, is of the KTB type. Here the gate 
charge has been found to play the role of an imaginary 
electric field, which grows under the renormalization- 
group transformation. In particular, it not only changes 
the nature of the transition but also tends to suppress 
quantum fluctuations, helping establish superconductiv- 
ity in the system. On the basis of this observation, we 
have obtained the zero-temperature phase diagram on 
the plane of the gate charge and the energy ratio, which 
is, unlike the diagram obtained in the previous analytical 
approaches, fully consistent with the results of quantum 
Monte Carlo simulations. 
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FIG. 1. Typical renormalization-group flow diagram on the 
K^ 1 — y plane. Here the initial (physical) value of the gate 
charge is q = 0.1 and the dotted line represents the locus of 
the initial configuration. 
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FIG. 2. Phase diagram of a Josephson-junction chain, dis- 
playing the boundary between the superconducting phase (S) 
and the insulating phase (I) on the Ej/Ec—q plane. The 
somewhat speculative boundary is depicted by the dashed 
line. 



